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Abstract The guided filter is a technique for edge-aware 
image filtering. Because of its nice visual quality, fast speed, 
and ease of implementation, the guided filter has witnessed 
various applications in real products, such as image edit¬ 
ing apps in phones and stereo reconstruction, and has been 
included in official MATLAB and OpenCV. In this note, we 
remind that the guided filter can be simply sped up from 
O(N) time to 0(N/s 2 ) time for a subsampling ratio s. In 
a variety of applications, this leads to a speedup of> 10 x 
with almost no visible degradation. We hope this accelera¬ 
tion will improve performance of current applications and 
further popularize this filter. Code is released. 

1. Introduction 

The guided filter [1, 2] is one of several popular algo¬ 
rithms for edge-preserving smoothing 1 . Its time complex¬ 
ity is O(N) in the number of pixels N, independent of 
the filter size. The guided filter can effectively suppress 
gradient-reversal artifacts [1] and produce visually pleas¬ 
ing edge profiles. Because of these and other properties, the 
guided filter has been included in official MATLAB 2014 2 
and OpenCV 3.0 3 and widely adopted in real products. 

Despite its popularity and its various third-party imple¬ 
mentations, we notice that a simple but significant speedup 
has not been exploited. This speedup strategy was briefly 
mentioned in [2] for joint upsampling but not for other 
generic scenarios. This method subsamples the filtering in¬ 
put image and the guidance image, computes the local lin¬ 
ear coefficients, and upsamples these coefficients. The up- 
sampled coefficients are adopted on the original guidance 
image to produce the output. This method reduces the time 
complexity from O(N) to 0(N/s 2 ) for a subsampling ratio 
s. An actual speedup of >10x can be observed. 

In this note, we revisit this speedup method by provid¬ 
ing more technical details, visual examples, and publicly 
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released code 4 . This acceleration method is particularly fa¬ 
vored for mega-pixel images, for which the filter size is usu¬ 
ally set as proportional to the image size in practice. As 
such, a local window on the subsampled images can still 
provide enough pixels for computing local statistics. In our 
extensive real applications for image processing, we have 
found that this speedup method has almost no visible degra¬ 
dation. Considering the growing usage of the guided filter 
in real products, we hope this simple speedup will improve 
the performance of these applications and further popular¬ 
ize this filtering technique. 

2. Method 

We denote the guidance image, filtering input image, and 
filtering output image as /, p and q respectively. The guided 
filter is driven by a local linear model: 

Qi = a k Ii + b k ,\/i G l c k , (1) 

where i is the index of a pixel, and k is the index of a local 
square window uo with a radius r. Given the filtering input 
image p , minimizing the reconstruction error [1] between p 
and q gives: 

R liPi - TkPk 

a k = —-2 i- ( 2 ) 

bk = Vk—^khk- ( 3 ) 

where p, k and a k are the mean and variance of I in the win¬ 
dow fc, and e is a regularization parameter controlling the 
degree of smoothness. The filtering output is computed by: 

qi = aili + bi , (4) 

where hi and bi are the average of a and b respectively on 
the window uji centered at i. The main computation is a 
series of box filters. Algorithm 1 shows the pseudo-code 
[2] of the guided filter, where / mea n( , ,0 denotes a mean 
filter with a radius r. 

4 http://research.microsoft.com/en-us/um/people/ 
kahe/eccvlO 
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Algorithm 1 Guided Filter. 

1 : mean/ — /mean (GO 
mearip - /meaner) 

COrr/ — /mean (7• * G 0 
COrr/p = /mean(/. * P, r) 

2 : var/ = corr/ — mean/. * mean/ 
cov/p = corr/p — mean/. * mean p 
3: a — cov/p ./ (var/ + e) 
b — mean p — a. * mean/ 

4: mean a — /mean (a, r) 

mean fe — /mean (GO 

5: q = mean a . * / + mean^ 


Algorithm 2 Fast Guided Filter. 

1 • I — /subsample (G £>) 

P — /subsample (p 7 s) 

r' — r/s 

2 : mean/ : — /mean(^ 7 r ) 
meanp - /meaner') 

COri/ = /mean(^ . * / ,7* ) 

COrr/p = /mean/ . * p ,T ) 

3: var/ = corr/ — mean/. * mean/ 
cov/p = corr/p — mean/. * meanp 
4: a = cov/p ./ (var/ + e) 
b — meanp — a. * mean/ 

5: mean a = / me an(a,r / ) 
meanj, — /mean (G r ) 

6 ; mean a — /upsample (mean a , s) 
rnean^ — /upsample (meant), s) 

1: q = mean a . * / + mean?, 


In the above, a and b in Eqn.(4) are two smoothed maps, 
and the edges and structures in q are mainly given by mod¬ 
ulating the image / (thus called guidance). But the major 
computation of guided filter is for the smoothed maps of a 
and b , which need not be performed in full-resolution. Al¬ 
gorithm 2 describes the subsampled version for Fast Guided 
Filter. We subsample (nearest-neighbor or bilinear) the in¬ 
put p and the guidance I by a ratio s. All the box filters are 
performed on the low-resolution maps, which are the major 
computation of the guided filter. The two coefficient maps a 
and b are bilinearly upsampled to the original size. Finally, 
the output q is still computed by q m a I + b. In this last 
step, the image / is the full-resolution guidance that is not 
downsampled, and it will still faithfully guide the output. 

The computation of all box filters reduces from O(N) 
complexity to 0(N/s 2 ). The last bilinear upsampling and 
output steps are O(N) complex, but only take a small frac¬ 
tion of overall computation. In practice, we have observed 
a speedup of >10x when s = 4 for both MATLAB and 
carefully optimized C++ implementation (depending on the 
number of channels). Fig. 1-4 show the visual results of 
using <s = 4 for various applications in [1]. 



image original GF Fast GF (s= 4) 

Figure 1. Edge-preserving smoothing, r — 4, e = 0.2 2 . The 
subsampling ratio is s = 4. 



original GF Fast GF (s= 4) 


Figure 2. Detail enhancement, r — 16, e = 0.1 2 . The subsam¬ 
pling ratio is s = 4. Top: edge-preserving smoothed images. Bot¬ 
tom: enhanced images with x 5 detail. 



guide input original GF Fast GF (5=4) 

Figure 3. Flash/no-flash denoising. r — 8, e = 0.02 2 . The sub¬ 
sampling ratio is s — 4. 



guide input original GF Fast GF (s=4) 

Figure 4. Guided feathering, r = 60, e = 0.001 2 . The subsam¬ 
pling ratio is s = 4. 
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